Quantum phase-slips in Josephson junction rings 
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We study quantum phase-slip (QPS) processes in a superconducting ring containing N Josephson 
junctions and threaded by an external static magnetic flux <J?s. In a such system, a QPS consists 
of a quantum tunneling event connecting two distinct classical states of the phases with different 
persistent currents [K. A. Matveev et al., Phys. Rev. Lett. 89, 096802 (2002)]. When the Josephson 
coupling energy Ej of the junctions is larger than the charging energy Ec = e 2 /2C where C is the 
junction capacitance, the quantum amplitude for the QPS process is exponentially small in the ratio 
Ej/Ec- At given magnetic flux each QPS can be described as the tunneling of the phase difference 
of a single junction of almost 2n, accompanied by a small harmonic displacement of the phase 
difference of the other N — 1 junctions. As a consequence the total QPS amplitude u T i US is a global 
property of the ring. Here we study the dependence of fring on the ring size N taking into account the 
effect of a finite capacitance Co to ground which leads to the appearance of low-frequency dispersive 
modes. Josephson and charging effects compete and lead to a non-monotonic dependence of the 
ring's critical current on N. For N — > oo, the system converges either towards a superconducting 
or an insulating state, depending on the ratio between the charging energy Eo = e 2 /2Co and the 
Josephson coupling energy Ej. 

PACS numbers: 74.50.+r,74.81.Fa,73.23.Ra,85.25.Cp 



I. INTRODUCTION 
A. One-dimensional Josephson junction chains 

One-dimensional Josephson junction chains (ID J J 
chains) have received considerable interest recently. 
Their use has been proposed for the realization of a qubit 
topologically protected against decoherence^— for the 
realization of a tunable parametric amplifier in narrow 
frequency ranges,^ for the realization of a fundamen- 
tal current standard in quantum metrologyj2ii£ and for 
designing controlled inductive electromagnetic environ- 
ments in quantum circuitry! 11 : 12 

Homogeneous JJ chains of infinite length have been 
studied theoretically in the pastJ^r— Such chains con- 
sist of superconducting islands, separated by Josephson 
tunnel junctions. In this paper we consider JJ chains 
arranged in a closed geometry, Fig. Q] The electrostatic 
interaction between the metallic islands is modeled by a 
neighboring capacitance C and by a local ground capac- 
itance Co, with Ec = e 2 /2C and E = e 2 /2Co the cor- 
responding charging energies. Each Josephson junction 
can sustain a maximum supercurrent Ij — 2eEj/h; this 
defines the Josephson coupling energy Ej. We consider 
rings 

Previous theoretical studies^ 3 . - — predicted a 
superconductor-insulator phase transition when the 
ratio between the Josephson energy Ej and the charac- 
teristic charging energy is reduced below a critical value. 
Bradley and Doniach studied this phase transition for 
infinite JJ chains and for the two extreme opposite 
cases when one of the two capacitances is vanishing 
(Co = or C = 0)»i2 Korshunov investigated the 




Figure 1: Schematic representation of a superconducting 
ring threaded by a magnetic flux <E>b and containing N iden- 
tical Josephson junctions with a mutual capacitance C and 
a local ground capacitance Co. Ej is the energy scale for 
the Josephson coupling. The variables (tp n ,Qn) are, respec- 
tively, the condensate phase and the excess charge of the nth 
superconducting island. 



general case for arbitrary ratio C/Coi 14 ' 15 He found 
that the critical value Ej ' of the Josephson energy 
at which the system undergoes the phase transition 
equals to Ej C) = E f(C/C a ) where f(x) is a smooth 
and regular function of order one. In particular, for 

(c) 

the case Co = 0, we have Eq = oo = Ej ' and the 
system is an insulator for any value of Ej in agreement 
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with the result of Bradley and DoniachJ^ Subsequently, 
experimental studies on the finite temperature behavior 
of the residual resistance in long one-dimensional chains 
of SQUIDs reported a phase transition when reducing 
the Josephson energyjiZr— 

The theoretical results reviewed so far were obtained 
in the thermodynamic limit TV — > oo. A first attempt 
to go beyond this limit was undertaken by Matveev, 
Glazman and Larkin 2 ^ who studied quantum phase-slip 
(QPS) processes in a superconducting nanoring contain- 
ing a large, but finite number of Josephson junctions. 
Here a QPS consists of a quantum tunneling between 
two distinct classical states of the phases with different 
persistent currents circulating in the ring at given mag- 
netic flux^ This is a collective process which can be de- 
scribed as the tunneling of the phase difference of a single 
junction by almost 2n, accompanied by a small harmonic 
displacement of the phase difference of the other TV — 1 
junctions (see Ref. [22| and the explanation in Sec. Ill Bl) . 
Quantum tunneling is possible due to the finite junction 
capacitance,which plays the role of inertia. 

Matveev et al. predicted a strong reduction of the 
maximum critical current sustained by the ring with in- 
creasing ring size TV due to QPS processes. Very re- 
cent experiments on flux-biased rings containing a few 
Josephson junctions 2 ^ indeed reported a remarkable sup- 
pression of the maximum supcrcurrent as Ej/En de- 
creases, in agreement with the findings of Ref. |22| . In 
these devices, the effects of the capacitance to ground 
could be neglected since the ring's circumference was 
much smaller than the screening length A of the system, 
given by A = -k^JC/Cq. However, it is expected that, for 
JJ rings of intermediate circumference TV > A, the effects 
of the capacitance to ground can no longer be ignored. 

In this paper, we study a JJ ring of finite circumference 
and threaded by an external magnetic flux $b, Fig. [T] 
We go beyond the previous work of Matveev et and 
we take into account the collective nature of a QPS as 
well as the ground capacitance Co for calculating the QPS 
amplitude. We show that the interplay between the finite 
value of the ratio Cq/C and finite size effects gives rise 
to a non-monotonic dependence of the low-energy prop- 
erties on TV. We systematically discuss this interplay as 
well as its consequences for the QPS amplitude for flux- 
biased rings with arbitrary number TV > 5. We focus on 
the limit where the Josephson coupling energy Ej dom- 
inates over the charging energies Ec,Eo, such that the 
amplitude for QPS to occur is exponentially small in the 
ratio Ej/Ec- This fact allows us to focus on the anal- 
ysis of a single QPS event. Once the QPS amplitude is 
known, one can calculate the ring's low-energy spectrum 
as a function of the external flux $ g and hence obtain the 
maximum supercurrcnt / max that the ring can sustain. 



B. Summary of the main results 

Our main results are summarized in Fig. [2] where we 
show the dependence of I max - scaled to the classical 
value I c i = irlj/N found in the absence of QPS processes 
- as a function of TV for two relevant situations: Co = 
and C = C/2. 

As we will discuss in detail below, a QPS can be de- 
scribed as an 27r(l — 1/TV) winding of the local phase- 
difference occurring on one of the junctions, accompa- 
nied by a simultaneous small (harmonic) adjustment of 
the phases of the other TV — I junctions^ The amplitude 
for the quantum tunneling of the phase-difference on one 
of the junctions of the ring between two different minima 
of the Josephson potential is given b y 22 ' 24 " — 

^ = ^{^E c y^J-f^\ . (i) 



a) b) 




Figure 2: Schematic behavior of the maximum supercurrent 
/max in 1D-JJ rings as a function of TV, scaled to the classical 
value 7 c i, for different values of the ratio Ej/Ec for (a) C'o — 
and (b) Co = C/2 (see also Fig. [S] and Fig. [TT] for details). 

Consider first the case Co = 0, Fig. HJa). The other 
TV — 1 junctions form a bath of dispcrsionlcss harmonic 
oscillators, all having the same (plasma) frequency uj v = 

1 /2 

(8EjEc) /h. In order to satisfy the constraint imposed 
by the flux threading the ring at all times during the 
QPS process, the phase differences for the TV — I other 
junctions perform a small shift ~ 1/N, on time scales 
1/uip comparable to the one needed to perform the actual 
phase hopping. This adjustment gives rise to finite-size 
corrections for intermediate ring circumferences TV to the 
amplitude given by Eq. (QJ leading to a QPS amplitude 
vq — > v(N). Since any junction can act as a QPS center, 
the total QPS amplitude for the ring is given by ^ r i n g = 
Nu(N). 

In the limit TV 3> 1, finite size effects vanish and v(N) 
converges to the constant v$ so that f r ing increases with 
the length and the maximum supercurrent / max vanishes 
exponentially^ see Fig. [UJa). The system becomes a 
perfect insulator at TV = oo, in agreement with Refs. 

MM 

On the other hand, we find that the interplay be- 
tween charging and Josephson effects in finite systems 
leads to an enhancement of the effective QPS ampli- 
tude v(N) with decreasing ring circumference TV, v(N) ~ 

1 /2 

vo cxp[(Ej/Ec) /TV], thus reducing the maximum su- 
percurrent / m ax- Consequently, the maximum supercur- 
rcnt / max shows non-monotonic behavior as a function of 
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the ring circumference TV for sufficiently large values of 
the ratio Ej/E C - 

When the capacitance to ground Co is restored, the 
TV — 1 harmonic junctions interact directly between them. 
This leads to the appearance of an ensemble of TV — 1 dis- 
persive electrodynamic modes at frequencies below u> p , 
similar to the ones found in a standard LC-transmission 
line^i The tunneling phase couples to these modes in 
much the same way as a quantum particle to a har- 
monic bath in the Caldeira-Leggett models In par- 
ticular, the low-frequency modes with linear dispersion 

1/2 

hiUk ~ {8EjEq) ' irk/N give rise to a finite friction for 
the QPS dynamics in the limit TV — oo. 

At finite TV, the coupling with the low-frequency modes 
strongly affects the QPS amplitude. Indeed, we find that 



1/2 



Dc- 



u(N) ~ v /N a for TV > 1 where a oc (Ej/E ) 
pending on the value of a, v T ing = Nv(N) ~ TV 1- " ei- 
ther tends to zero, when a > 1, or grows linearly, when 
a <C 1, indicating that the system either displays a su- 
perconducting or an insulating behavior, as can be seen 
in Fig. Hfb). This behavior is reminiscent of the dissipa- 
tive phase transition^ occurring in a single junction in an 
electromagnetic environment ! 14 ' 15 For intermediate ring 
sizes finite size effects occur, yielding a non-monotonic 
behavior of the maximum supercurrent in the insulating 
regime, similar to what we find for the case Cq = 0. 

The paper is structured as follows. In Sec. [TT] we re- 
call the model for a flux-biased ID JJ ring as well as 
the notion of QPS and the approximations involved to 
find the ring's flux-dependent quantum ground state and 
hence the maximum supercurrent. In Sec. Mil we discuss 
the single QPS approximation and we show how the sys- 
tem reduces to a model similar to that of Caldcira and 
Leggett^ where one single Josephson junction, the cen- 
ter of the QPS, is coupled to TV — 1 harmonic oscillators. 
The results for the specific case Co = are shown in 
Sec. lIVI where we explain in detail the different finite-size 
corrections on the QPS amplitude. The effect of the fi- 
nite ground capacitance Co > and the general results 
are discussed in Sec. [V] In the last Section fVTl we draw 
our conclusions. 



II. THE MODEL 

A. Hamiltonian 

We consider a homogeneous ring of TV identical super- 
conducting islands each coupled to its nearest neighbor 
by Josephson tunnel junctions, sec Fig. [I] We assume the 
superconducting gap A of the islands to be the largest 
energy scale involved in the problem. If A 3> 6E where 
5E is the average spacing of the electronic energy lev- 
els, superconductivity is well established. The islands 
should be metallic with large enough volume so that the 
perturbative treatment of Cooper pair tunneling through 
the contacting surfaces is justified. We assume the ab- 
sence of quasi-particlc excitations at low temperature T 



and low voltage fceT, 2eV <C A, where V is the typical 
voltage across the junctions. Furthermore, the condi- 
tion A ^> Ec , Eq implies that the Josephson coupling 
energy Ej characterizing Cooper-pair-tunneling between 
islands is independent of Ec,Eq. These conditions uni- 
vocally define the standard quantum-phase model for 1D- 
JJ chains.— Moreover, although various boundary condi- 
tions are in principle possible, we will henceforth consider 
a flux-biased ring threaded by a magnetic flux &b (see 
Fig.[TJ. The Hamiltonian reads 



H = - ^ QnC n ^Q,n - Ej ^ COS (0 n+1 - n ) 



N-l 



N-l 



n,m— 
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kin 



2e 



N-l 
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,n=0 



n=0 



2n m 



$11 



(2) 



where $o is the quantum flux, m is an integer and the 
periodic index 6 n =N = 6n=o- For each island, the BCS 
condensate phase (p n and the excess charge Q n on the 
nth island represent the conjugate variables of the system 
[<fn,Qn] = 2e«. The relative phase-difference across the 
nth junction between neighboring islands n and n + 1 is 
8 n = (fin+i — ifin- C is the capacitance matrix with matrix 
elements C„, m = (C + 2C)<5„, m - C(<5„ + i, m + £„_i, m ), 
with the index n = — 1 corresponding to TV — 1. 

We assume that the kinetic inductance, associated 
with the kinetic energy of the Cooper pairs in each su- 
perconducting island, is negligible as compared to the 
Josephson inductanc o 30 ' 31 L/, in — s- so that the last term 
in Eq. © is equivalent to the constraint^ 



N-l 

71 = 



27rm + 2-7T 



<Fb 

$0 



(3) 



We also assume that the geometric inductance can be 
neglected so that the current circulating in the ring docs 
not generate any magnetic field. The total flux $s is 
thus only given by the externally applied magnetic field. 

From Eqs. © and we see that the physical prop- 
erties of the system depend periodically on the ratio 
S = 2tt<I>_b/$o- In the steady-state regime, the dc- 
supercurrcnt flowing through the ring is the same for all 
the junctions (/„) = I and can be related to the deriva- 
tive of the ground state energy Eos of the system with 
respect to 6 



1(6) 



dE t 



as 



9$, 



dE, 



GS 



(4) 



It is in general a difficult task to find the ground- 
state energy Egs(S) for the flux-biased ring described 
by Hamiltonian ([2J and ([3]). An approximate solu- 
tion can be found in the limit where the Josephson en- 
ergy is larger than the characteristic electrostatic energy 
Ej 3> Ec,Eq, as we will discuss below. 
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B. Single QPS in JJ rings of finite circumference 

To set the stage, let us first consider the classical limit, 
achieved by setting Eq = E = 0, so that the phases are 
well-defined classical variables. The classical energy of 
the system reduces to 



N-l 

Ed = ~Ej C0S ( 6 n) 
n=0 



(5) 



The energy Eq. ([5]) is invariant under a change by 2n 
of the phases 9 n (or <p n ). In other words, the state 9 n 
and 9 n + 2itm represent the same state. However, at 
fixed magnetic flux 6, a given configuration of {9 n } cor- 
responds to a real physical state only if the constraint 
Eq. is satisfied. Therefore any distribution of the 
phases that violates Eq. (J3J is unphysical. 

The classical states \m) that minimize the energy 
Eq. ([5]) under the constraint ([3]) correspond to a uni- 
form distribution of phase differences 9 n = (2irm + S)/N. 
They have energies 



E,, 



-E T N cos 



277777, + 5 

TV 



(0) 



with the condition d 2 E c i/d9^\ = —E m > and the 
index -(N-l)/2 < m < (N-T)/2 (N odd) or -N/2 < 
to < (N/2) — 1 (N even). These classical states |m) are 
physically distinguishable as they are characterized by 
different persistent currents 



/,/ N sin 



2irm + 5 
N 



(7) 



Away from the degeneracy points 8 = and 5 = tt, they 
have also different energies. The classical ground state 
corresponds to an absolute minimum of potential energy 



. 2nm + S\ Ej . 2 

Efiq = —EjN max cos ~ mm (2-Km + o) , 

m \ N J 2N m 

(8) 

where the second, approximate equality is numerically 
accurate for sufficiently long ring (N > 5) so that the 
corresponding supercurrent has as a sawtooth-like depen- 
dence as a function of S with a maximum supercurrent 
given by I cl ~ %Ij/NM 

For finite C, Co , the electrostatic interaction acts as an 
inertial term on the phases so that quantum fluctuations 
occur giving rise to quantum phase-slip (QPS). At fixed 
magnetic flux and in a ring of finite circumference, the 
QPS is a collective process corresponding to the quantum 
tunneling in a multidimensional space of dimension N be- 
tween two distinct minima of the potential, correspond- 
ing for instance to the classical states \m) and \m + 1), 
separated by some energy barrier associated with the po- 
tential Eq. ([5]). In the multidimensional space, the physi- 
cal paths {9 n } which connect \m) and |m ± 1) correspond 
to a subspacc defined by the constraint Eq. J3J). Due to 



this constraint, the multidimensional tunneling reduces 
to one-dimensional tunneling in which we have only a 
few trajectories connecting the initial state and the final 
state, see Fig. [3J As it was discussed in Ref. [H, an ex- 
ample of QPS connecting the states \m) and \m + 1) is 
given by the displacements 



5 + 2nm 
N 

5 + 27rm 
N 



<5 + 27T(?77- 1) 

N 

5 + 2ir(m- 1) 
N 



f 2tt (9) 
(77 ? n ) , (10) 



in which the local phase difference 9 no , the center of the 
QPS, winds by an amount of 27r(l — 1/N) and the whole 
set of phase differences {9 n } (n ^ no) shifts in order to 
preserve the constraint. Fig. [3] shows this process. One 
can express the classical energy Eq. ([5]) as E c i/Ej = 
- cos(0„ o ) - (N - 1) cos(0 n ) with 77 ^ n , Fig. [3p. 

As it is shown in Fig. [3Jd, the trajectory of the QPS 
Eqs. (|9"j) . (fTU|) can be drawn in a restricted region as well as 
in an extended one. This fact allows us to introduce the 



a) 




Figure 3: Example of a QPS process, a) The oval loop 
represents the JJ ring. The arrows represent a few phase- 
differences {6 n } around the QPS center, the phase 8 no wind- 
ing of 2n(l — 1/N) (red arrow). The initial configuration 
(left) is the state \m) and the final configuration (right) is 
]m + l). The central configuration is intermediate between 
the two states \m) and |m + l). b) The classical energy 
expressed as E c i/Ej = — cos(x) — (N — 1) cos(y) with the 
axis x = 9 na and y = 9„(n 7^ no). The dots represent the 
state \m) and |m + l). The cross also represents the state 
\m + 1) but in the extended region. The black lines represent 
the physical trajectory of the QPS that connects the initial 
and the final state, the (black) solid one in the restricted re- 
gion, the (red) dashed line in the extended region. The bold 
line for the borders marks the compact region. Parameters: 
AT = 10,5 = 7r/2, m = 1. 
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adiabatic potential for the QPS process. In the limit in 
which the evolution of the phases is extremely slow, the 
kinetic energy can be neglected at any time and one ob- 
tains the energy of the system just minimizing the classi- 
cal energy for each intermediate configuration which con- 
nects the initial and the final state. Hence the adiabatic 
potential is associated with the line of minimum energy 
on the surface Eq. ([5]) which connects the end points as 
shown in the example of Fig. [3Jd. This line is given by 
the condition 



2irm — 9 r , 
TV- 1 



(n ^ n ) , 



(11) 



so that the adiabatic potential V e ff(0 no ) reads 

S + 2inn — 9 r - 



V eff(O no ) = -Ej 
~ Ej 



cos(0„ o ) + (TV-l)cos 

(5 + 2irm-0 no ) 



TV - 1 

2" 



-cos(6»„ ) + - 



2(TV- 1) 



(12) 



where the second, approximate equality is valid for suf- 
ficiently long rings (TV > 5). At the initial time we have 
V e ff [9 no = (5 + 27rm)/TV] = E. m and at the final time 
Veff{O no =2tt+[5 + 2v{m + 1)]/N} = E m+1 . 

V e f f (9 no ) is shown in Fig. 2] In the limit of TV — > oo 
we have V e ff(8 no ) = — Ej cos(#„ ) and we recover the 
simplified picture of QPS corresponding to the quantum 
tunneling of the phase difference 9 na of a given junction 
no from one minimum of the local Josephson potential 
— cos(#„ n ) to the neighboring one. 




1 s~ \ x -2it 271 4n J i i 

\\\ /<NA 

V / x ^ ' 8=71/2 



Figure 4: Example of the effective potential V e ff{0n o ) 
Eq. |12J for TV = 10 and for the values 5 = solid (red) 
line, S — 7r/2 dashed (purple) line and 5 = n dotted 
(blue) line. Inset: comparison between the effective poten- 
tial at <5 = 7T, dotted (blue) line, with the cosine potential 
EjCOs[0 no (N - 1)/N], solid line (see text). 

On the other hand, as we explained above, QPS is a 
collective process corresponding to a quantum tunneling 
in a multidimensional space with the constraint Eq. ([3]) 
imposed by the magnetic flux threading the JJ ring. Con- 
sequently, the potential in terms of the variable 9 no , the 
center of the QPS, is not a periodic function although 



the global classical energy of the system Eq. ([5|) is 2ir 
periodically invariant. 



C. Effective low-energy description 

We now turn to the effects of quantum phase fluctu- 
ations in the limit Ej 3> Ec,Eq. In this case, a simple 
analysis is possible since the QPS processes occur only 
rarely, with an exponentially small amplitude v, Eq. ([1}. 

As a consequence of QPSs and for small amplitude 
v, the quantum ground state I^gs) of the ring corre- 
sponds to a superposition of different classical states of 
the phases, namely I^gs) = H m c m |m). The coefficients 
Cm as well as the quantum ground state energy Egs (8) 
can be obtained from the following effective Schrddinger 
equation^ 



E n 



'ring 



-i) — Eqs c n 



(13) 



where the term proportional to Z4-i ng = Nv connects 
two classical states differing by a single QPS. The fac- 
tor TV takes into account the fact that the center of a 
QPS can occur at any junction corresponding to differ- 
ent trajectories in the multidimensional space TV. In 
other words, QPS amplitudes add up coherently. This 
coherence has been recently confirmed experimentally in 
a short 6-SQUID JJ chain^ by the measurements of the 
Aharonov-Casher interference effect. The coherence is 
affected by off-set charge dynamics. The details of this 
dynamics are currently not understood. It is expected to 
give rise to an additional dependence of v on TV which is 
beyond the scope of the present paper. 

The behavior of the general solution for the ground- 
state of the model given by Eq. (fl~3"|) is determined by 
only one dimensionless parameter— 



TV 2 ^(TV) 
2tt 2 E t 



(14) 



The ground-state energy Egs = Egs '(8; q) depends para- 
metrically on q, hence leading to a dependence of the 
maximum supercurrent I ma ^{q) on the QPS amplitude 
v. This dependence is illustrated in Fig. [5] together with 
the evolution of the value of the phase S max that cor- 
responds to the maximum. For q <C 1, the solution 
scales as I max /Id — 1 — (5/2)g 2 / 3 whereas it scales as 
I m ax/Ici ^ 24 v ^Fg 4 / 3 exp(-8 v ^) for q > 1. 

The effective low-energy theory defined by (|13|) was in- 
troduced in Ref. [13 for Co = and was discussed for long 
J J rings TV 3> 1. As we discussed in sec. Ill B[ the valid- 
ity of Eq. (|13|) extends more generally to include rings 
of intermediate size and with Co ^ 0. In the follow- 
ing, we will obtain the detailed dependence of the QPS 
amplitude v T - ms on the parameters TV, Ej,Ec, and Eq. 
We calculate the ground state energy of the system using 
Eq. (jT3")) . This procedure allows us to obtain the periodic 
dependence of Egs and I ma x on 8, Eq. (HJ). 
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0.2 0.4 0.6 0.8 

q 

Figure 5: Maximum supercurrent (solid line) for the model 
defined by Eq. (fl3)l as a function of q, Eq. (fl4|) . The dashed 
line represents the value of the scaled flux <5 max which gives 
the maximum supercurrent. 



where the Euclidean action for the phases {<^n(r)} reads 
S = J drC, and L is the Lagrangian, 



n=0 n=0 



N-l , g 

^2 E J COS ( (f n+ i ~ ( Pn + ^ 

71 = ^ 



(16) 



with tp n = dtp/dr and 5 m = S + 27rm. The compact vari- 
ables {<^n} are defined on the circle [0;27r[. Notice that 
we shifted the phase differences {9 n } with respect to their 

average value so that Eq. ([3]) now reads J2 n =o ® n = 0- 
The last constraint is automatically satisfied by imposing 
the boundary condition ip^ = ipQ. 



III. QPS IN ID JJ RING 
A. General approach for a single QPS event 

We present a general approach to calculate the quan- 
tum amplitude v for a single QPS event occurring on a 
ring containing N junctions and with mutual and ground 
capacitance C and Co, in the regime Ej ^> Ec, Eq. 

Recalling Eqs. (|9I10[) . the QPS process is a collective 
process in which the local phase difference 8 no winds by 
an amount of 27r(l — 1/N), accompanied by a shift of 
the whole set of phase differences {0 n } (n ^ no). For 
rings of circumference N > 5, the phase differences of 
the other junctions will vary only slightly A8 n ~ 1/N 
compared to the period 27r of the cosine potential. Hence 
we can apply the harmonic approximation to describe the 
dynamics corresponding of 9 n , n =/= uq . 

For vanishing ground capacitance Co = 0, the phase- 
differences of the other N — 1 junctions behave as in- 
dependent LC— oscillators at the plasma frequency lo p 
whose displacement is inversely proportional to the cir- 
cumference N. Thus, in the limit of very long ring 
N 3> 1, the dynamics of the other N—l phase differences 
can be neglected^ so that the QPS amplitude v(N) can 
be approximated by the N— independent constant value 
isq , Eq. ((T|). However, for ID- J J rings of finite circum- 
ference, the dynamics of the other junctions can have 
considerable effects on v, as we will show below. More- 
over, for finite ground capacitance Co > the harmonic 
oscillations of the phase-differences of the TV — 1 junctions 
play a crucial role for any ring's circumference. 

To calculate the QPS amplitude v, we start by consid- 
ering the partition function associated with the Hamil- 
tonian Eq. ([5]), with Lki n — > giving the constraint, 
Eq. ((3j). In the path integral formalism, it reads (/? = 
h/k B T) 



N-l 



Tr 



-H 



71 = 



[ fV[<p n (T)] 



-S/h 



(15) 



B. Harmonic modes 

Let us briefly discuss the behavior of the system in 
the harmonic approximation, neglecting the QPS. When 
the Josephson energy Ej ^> Ec,E the phases fluctu- 
ate only slightly around their classical values. The aver- 
age phase difference between the neighboring islands is 
small so that we can expand the Josephson interaction 
to lowest (quadratic) order. The general imaginary-time 
Lagrangian Eq. (|16|) then reduces to the harmonic one, 
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where we omitted an irrelevant constant term. Any 
periodic function <p n defined on the finite lattice n = 
0, . . . , N — 1 can be decomposed as 
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with the condition for the complex variables tfN—k = 
which preserves the total number of degrees of freedom. 
Substituting Eq. (Q2J into Eq. (JTTJ) and summing over the 
index n, the harmonic Lagrangian is diagonalized (see 
also Appendix [A"|l 
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the frequency dispersion is given by 
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where the screening length is A = ir^fC]C~ a . For N > 1, 
the maximal frequency of the modes uj max is given by 
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An example of this frequency dispersion is given in Fig. El 
We can distinguish two regimes for Cq > 0. For JJ rings 
longer than the screening length A, N 3> A, the spec- 
trum has a linear dispersion for low frequencies. For 
shorter rings N <C A, the lowest mode appears almost 
at the plasma frequency uj p and the linear behavior is 
completely absent. For the case Co = we recover a 
flat distribution where all the modes are degenerate and 
correspond to the plasma frequency uj p . 




Figure 6: Frequency dispersion of the harmonic modes in a 
JJ ring of N — 100 junctions. We use the equivalent index 
notation k = —N/2 + 1, . . . , N/2. Solid (purple) line, dashed 
(blue) line and dotted (red) line are, respectively, for the ratio 
C/Co = 0.1,1,100. 

Let us calculate the phase-dependent ground state en- 
ergy Eqs* ^ (5) in the harmonic approximation. The path 
integral in Eq. (|15[) can be explicitly calculated for the 
diagonalized harmonic action Eq. (IT9l) . We obtain 
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where the product corresponds to the partition function 
for an ensemble of N— 1 independent harmonic oscillators 
and the exponential factor contains the classical energy 
of the system (mod 27r) at finite temperature. Notice 
that the zero frequency mode k = is not involved in 
the relevant part of the partition function. Note also 
the periodicity of the result Eq. (|23|) with respect to S. 
The quadratic dependence on 6 in each of the segment 
— 7r + 2irm < S < —it + 2ir(m + 1) is a consequence of the 
harmonic approximation, i.e. Eq. ([8]). 

We thus obtain the phase-dependent ground-state en- 
-lim p=00 (h/(3)\nZ^ 
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Notice that to reach the full quantum regime of the har- 
monic modes, as expressed by Eq. (f2"4"|). the temperature 
has to be much smaller than the energy of the lowest 
mode, ksT <C which represents the lowest en- 

ergy scale of the system given by 
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We conclude that, in the absence of QPS processes, the 
JJ ring forms a closed quantum LC-line (each Joseph- 
son junction in Fig. [T] is replaced by an inductance 
Lj = ft 2 /Ae 2 Ej) formed by N — 1 independent harmonic 
oscillators with eigenfrequencies u>k- The classical saw- 
tooth relation between the supercurrent and the total 
phase I = 1(5) is unmodified by the quantum harmonic 
fluctuations of the phases. However, harmonic fluctua- 
tions are relevant when the QPSs are restored, as we will 
now show. 



Effective QPS action in the presence of 
harmonic modes 



As discussed in Sec. Ill Bl when a single QPS is centered 
in one specific junction no, the dynamics of the other 
junctions (n ^ hq) is well described by the harmonic ap- 
proximation. Then, as seen in Sec. lIIIBl we can consider 
this part of the ring as an electromagnetic environment 
formed by N — 1 independent harmonic oscillators and 
to which the single junction, center of the QPS process, 
is coupled. Accordingly, it is natural to cast the full ac- 
tion in a form where the winding phasc-diffcrcncc 9 no is 
linearly coupled to an ensemble of harmonic oscillators 
acting as an external bath. 

Without loss of generality, one can assume the center 
of the QPS to the junction n = N — 1, namely 8 = 
9n-i = <^o — <Pn—i- It is useful to introduce the average 
value of the phase at the junction N — 1, namely <po = 
((po + tfN _i)/2, so that one can write 
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We choose as the relevant variables the set S given by the 
winding phase difference difference 8 together with the 
local phases, the average tpo and (p n with n = 1, . . . , N—2. 

First we discuss the harmonic expansion of the poten- 
tial energy of the Lagrangian Eq. (|TB)) . In a QPS, the 
phase-differences across the junctions remain small with 
respect to 2tt except at the junction n = N — 1 
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Using the set S, the last sum in Eq. (|27|) contains a 
quadratic coupling between 8 and the quantities if>Q — ifi 
and tpo — tfN-2- Although these terms contain 8 varying 
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by almost 27r, this does not make invalid our expansion as 
the overall argument of the cosine function, representing 
the phase-difference in the neighboring junctions n = 
and n = N — 2, remains small during the QPS. 

Using the set S, the Lagrangian of Eq. (fT6|) is decom- 
posed as 



C = Ci + C 2 + C 3 
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where the first term C\ is associated with the winding 
junction 
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The second term £2 describes the environment to which 
9 is coupled and it reads 
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with the notation ip a = ip^-i = <fio for the periodic 
boundary condition. Note the extra term in Eq. (|3"0")) 
associated with the phase at n = 0. The last term £ 3 
of Eq. ([28]) describes the coupling between the winding 
junction and the electromagnetic environment, 
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We rewrite the Lagrangians £2 and £3 in terms of the 
normal modes tp k . They are given by Eq. ([TB"]) with N 
replaced by N — 1. In terms of these modes, we have 
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where X) is defined by Eqs. (19]), ([20]), dHJ in which 
we have to replace iV by N— 1. Then, from Eq. (f3"Tj) . one 
can see that the phase 9 is coupled only to the imaginary 
part of the modes ipk = <p k + ( see Appendix \K§ 
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where k max = (N — 3)/2 for A^ odd and k max = 
for N even. We have introduced the factor 

C fe = 4 sin [27rfc/(iV - 1)] /(TV - 1) 1/2 
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Notice that the mode fc = is not coupled to the wind- 
ing phase 9 and we can omit it hereafter. The partition 
function associated with the total Lagrangian Eqs. ([28]) , 
([29]). ([32]) . ([S3]) reads 
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It is possible to integrate out the imaginary parts of the 
N — 2 harmonic modes to obtain a single effective ac- 
tion describing the dynamics of 9. After the calculation, 
shown in Appendix fX] we find 
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where Zj^ r 2 ^ is given by Eq. ([23]) with A^ — 1 replacing 
iV. After a shift of the phase 9 — » 9 — 6 m /N, the effective 
action for the phase is given by 
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The effective action has a kernel G(r) which is non-local 
in time and whose Fourier series is given by G(t) = 
(l/P)Gg exp(iw^r) where = 2n£/f3 are bosonic 
Matsubara frequencies and G^ reads 
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The kernel has the relevant property Gi — for I = 
W£ = 0. The last relation is equivalent to J^(1tG(t — 

t') = J^cIt'G(t — t') = 0. As a consequence, this 
term is invariant under a shift of the winding phase 
9 — > 9 + const. In other words, upon a proper redefinition 
of G(w) = lu 2 G'(u>), this (kinetic) term can be written as 
~G'(t-t')9(t)9\t'). 

We observe that the potential in the first line of S e ff 
in Eq. (|37[1 corresponds exactly to the adiabatic potential 
V eff {9) Eq. (TTJI introduced in Sec. iHTfl This potential 
breaks formally the periodicity in 9 in the action. This 
symmetry breaking is a consequence of the fact that the 
QPS is a quantum tunneling in a multidimensional space 
with the constraint imposed by the magnetic flux thread- 
ing the J J ring (see discussion in Sec. IIIBI) . 

In summary, Eqs. . ([5T]) . and ([33]) constitute the 
central result of this paper. They enable us to calculate 
the size-dependent QPS amplitude v(N) and hence the 
phase-dependent ground-state energy and the ring's su- 
percurrent in a broad range of values of the parameters 
N, Ej, Ec and Eq, as we will show in detail below. 

However, we first establish a relation with previous 
wor k 1 15 on infinitely long chains by considering the 
thermodynamic limit. 



D. The thermodynamic limit and the dissipative 
dynamics 

The effective action Eq. ([3T]) describes the single wind- 
ing junction coupled to its electromagnetic environment 
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constituted by the other N — 2 junctions in the harmonic 
approximation. This action is very similar to the one 
describing the dissipativc dynamics of the single Joseph- 
son junction in the framework of the Caldeira-Leggett 
models In this model an abstract bath formed by an 
infinite number of harmonic oscillators is phenomenolog- 
ically introduced as the mechanism of irreversible loss of 
energy in the Josephson junction. 

The external bath discussed here, expressed by the 
kernel G(r) in Eq. ([55]). physically corresponds to the 
reai harmonic modes sustained by the Josephson junc- 
tion ring. These discrete modes can be experimentally 
designed and tested*^ As long as the ring has finite size, 
there are a finite number of discrete modes and no real 
dissipation occurs. 

We also note that the interaction between the wind- 
ing local phase-difference at the junction no and these 
N — 2 harmonic modes is characterized by a linear cou- 
pling through the positions of the oscillators (fk as well 
as through their velocities (fk, see Eq. 

As we will show now, the difference between the system 
described by Eqs. (|5Tj) . and the standard Caldeira- 
Leggett model disappears in the limit N = oo. 

Let us consider Eqs. ([3T [> .(|38 ]) . Taking the limit N -> 
oo, the first term on the right hand side of Eq. (|5Tj) 
reduces to the action for a single capacitively shunted 
Josephson junction with a capacitance C + Co/2. Equa- 
tion (|38|) for the kernel Gi then takes a simple form by 
replacing the sum with an integral. Proceeding in this 
way, we add the kinetic term of the first line of Eq. ([57)1 
to G( to recover Korshunov's result for the total kernel 
7(0;^ ) of a QPS in a chain of infinite length^ 

K^iaxJ 4 V 2£ V ^fnax 

where u) max is defined in Eq. (f22j) . The effective action 
now reads: 

rP rP i rP 
S e ff= dr /dT , -7(T-r')0(T)0(T')- dTEjcos[0{T)] . 
Jo Jo * Jo 

(40) 

First, let us discuss the high and low energy regions 
for the dynamics of the winding phase difference 9 in 
the cosine Josephson potential. These two regions are 
separated by the condition that the kinetic energy be re- 
spectively larger or smaller than the height of the poten- 
tial ~ Ej. To estimate the kinetic energy, we determine 
the effective capacitance of the junction. This can be 
achieved by taking the limit ujg. — * oo. Then the kernel 
j(u>i) Eq. (|39j) corresponds simply to a pure capacitance 




The kinetic energy corresponds to the height of the po- 
tential barrier when 7(0^) = Ej giving the threshold 



uj£ = uimax- For cj£ > ui m ax the cosine Josephson poten- 
tial is a small perturbation for the dynamics 9. On the 
contrary, the range uj£ < U3 max corresponds to the energy 
region where the winding phase 9 moves well within a po- 
tential minimum. In the absence of the interaction with 
the modes, its dynamics is harmonic. 

The threshold uj max separates the low and high- 
frequency part of the kernel 7 Eq. (|39|) . As we have 
explained, the high-frequency component u>£ > u> m ax 
scales approximately as 7(w<) ~ and the chains be- 
haves as a pure capacitance. On the other hand, ap- 
proaching the zero frequency, 7(0^) scales linearly. In a 
Josephson junction with dissipation, this behavior cor- 
responds to the effect of an external resistance R lead- 
ing to a dimensionless damping parameter 77 = R q / R 3 ° 

where R q = ft/4e 2 . In our model, for N — 00, this resis- 

1/2 

tance corresponds to R — (Lj/Cq) , the characteristic 
low-frequency impedance of the chain related to linear 
dispersion of the modes. 

It is interesting to observe that in the standard 
Caldeira-Leggett model the two energy scales -the first 
one related to the ratio between the height of the poten- 
tial and the kinetic energy and the second one related 
to the high-energy cut-off for the dissipation- are gen- 
erally different as the latter is determined by the spe- 
cific electromagnetic environment to which the junction 
is coupled. 

Defining the relaxation time r r = RC, we can calculate 
the quality factor given by the ratio between the oscil- 
lation period and the relaxation time, Q = 2ir/(ujpT r ) = 

1 /2 

2h(C/Cq) . In Josephson junction chains, the under- 
damped regime Q ^> 1 corresponds thus to C ^> Cq. 



IV. RESULTS: VANISHING GROUND 
CAPACITANCE 

We now go beyond the thermodynamic limit and dis- 
cuss finite-size effects for a ring in the limit Co = 0. The 
non-local kernel then vanishes, Gi = 0, and the effective 
Lagrangian from Eq. ([57)) reduces to 

= \— ) 8^ " Ej + 2{N-1) ■ 

(42) 

This result has a simple interpretation. Let us write the 
total Lagrangian, Eq. ([28)) . with Co = in terms of the 
phase differences 9 n , 

We observe that the phases {9 n } are not coupled ex- 
plicitly to the tunneling phase 9. But the total phase 
difference is fixed 9 n + 9 = 5 m which enforces an 
implicit interaction between them. In case of identical 
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phase-differences 8 n for all N — 1 junctions, we have 



of the energy barrier separating the two wells), 



Sjn - 

N - 1 
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(44) 



By inserting Eq. (|4"4")l in the action Eq. (|4"5|). we obtain 
the effective action Eq. (|4"2j) . Clearly, in the limit A" — > oo 
finite size corrections vanish and we recover the simple 
action for the single QPS. 

A finite value for N first of all modifies the kinetic 
term. The constraint-induced coupling to the other junc- 
tions increases the inertial mass C* jC of the single phase 
performing an almost complete winding, leading to a re- 
duction of the charging energy, 
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(45) 



This causes a reduction of the QPS amplitude. 

On the other hand, as we can see from Eq. ([4"2]). the 
action of the phase 8 involves the effective adiabatic po- 
tential V eff (9)/Ej = l-cos8 + (6-6) 2 /{N- l)/2, see 
Fig. which is not purely sinusoidal and depends on 
N, leading to additional finite-size corrections. 

Before discussing these corrections, we note that for 
short rings, the effective potential can have only one 
minimum for arbitrary value of 6. This is the case, for 
instance, when S = for N < 5. In this regime the 
harmonic approximation fails too. Focusing on lengths 
longer than N > 5, the potential has always two minima 
for which the QPS process is well defined. 

Due to the periodicity in S m and the symmetry in 
(5 m <H- — 5 m ) for the quantum ground-state of the sys- 
tem, we can restrict our discussion to values of S be- 
tween and 7r (m = 0). In this case, the relevant tun- 
neling process corresponds to an increase of the variable 
8. The two positions corresponding to the two minima 
of the effective potential are, respectively, 8i = 5/N and 
8 r = 2tt + (5 — 2k) /N, with the energies V e ff(0i) and 
V e ff(8 r ) which correspond to the two classical energies 
associated with the initial state and the final state before 
and after a QPS, respectively. 
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The ratio Ej/Ej is an increasing function of N which 
converges to one in the limit of N ^> 1. Thus for any 
finite-size system, this correction corresponds to a de- 
crease of the barrier for quantum tunneling which there- 
fore would lead to an enhancement of the QPS amplitude, 
in contrast to the effect of the renormalized capacitance, 
Eq. (|4"5|) . Another effect which enhances the tunneling 
amplitude is the reduction of the distance between the 
two minima of the effective potential, 



A8* _ 1 
~2~k ~ ~ ~N ' 



(47) 



contributing to the enhancement of the tunneling ampli- 
tude for finite N. Similar renormalization effects have 
been discussed in Ref. H^, albeit for different supercon- 
ducting circuits. 

The three different effects given by Eqs. (|45l) . (J46J , 
(|47p combine in the final expression for the renormalized 
effective amplitude for the quantum tunneling between 
the two degenerate minima, 
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(48) 



For the range of interest N > 5, the QPS amplitude 
Eq. (|48|) decreases with the length N indicating that 
the latter effects (reduced barrier height and distance 
between minima) dominate the capacitance renormal- 
ization. Indeed, when v is scaled with the amplitude 
Vq (corresponding to the limit A^ = oo), the leading 
exponential term reads v/vq ~ cxp(N c /N) with N c = 
(it 2 + 4) yJ%E,jjE C - Here A^ is the typical length below 
which this finite size correction becomes relevant. The 
behavior of the ratio v/vq is shown in Fig. [7] 

The effective QPS amplitude v decreases with increas- 
ing N. As a result the parameter q ~ N 2 v(N) defined 



A. Symmetric effective potential (q < 1) 

For small values of the parameter q Eq. (|14p , the hop- 
ping term is small, Eq. (|13|) . The QPS process couples 
mainly two neighboring classical states, |m) and \m + 1), 
and it is significant when they are degenerate. This oc- 
curs for S ~ 7T which corresponds to the point of the 
maximum of the supercurrent in this regime (Fig. [5]). 

Therefore, we start by discussing the case S = ir when 
the effective potential V e ff(8) is symmetric around the 
maximum, Fig. ((4]). Then the effective potential is very 
well approximated by a renormalized cosine potential of 
the form -E} cos{8(N- 1) /N] (see inset Fig.HJ) where E} 
is the renormalized Joscphson energy (half of the height 
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Figure 7: For Co = 0, the renormalized QPS amplitude v as 
a function of N, Eq. (|48[l . scaled with i>o, Eq. ^Q, for different 
values of the ratio Ej/Ec- 
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in Eq. (JT1} has 

of N. According 
for Eq. ((T3J) has 
total length N. 
percurrent as a 
ratios of Ej /Ec 
behavior occurs 



a non-monotonic behavior as a function 
dy, the maximum supercurrent obtained 
also a non-monotonic dependence on the 
In Fig. [5] we show the maximum su- 
function of the ring size N for different 
Decreasing Ej j Ec, the non-monotonic 
at shorter lengths. 
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Figure 8: The maximum supercurrent I ma x scaled with / c i as 
a function of N and at different values of the ratio Ej /Ec by 
using the QPS amplitude u(N) Eq. (|48[) . Inset: the maximum 
supercurrent by using the QPS amplitude vq, Eq. {TJ. 



B. Asymmetric effective potential (q > 1) 

When the parameter q is of order ~ 1, the maxi- 
mum supercurrent for the model given by Eq. (|13[) shifts 
from the phase difference 6 = ir to the phase differ- 
ence S = 7r/2, Fig. [5J The current-phase relation is 
strongly modified passing from a sawtooth to a sinusoidal 
function^ 

We should then solve Eq. (fP3")) using the renormalizcd 
QPS amplitude v as associated with the full asymmetric 
effective potential V e /f and which depends on S, Fig.|4j 

To reach the regime <? > 1 there are two possible ways. 
First, it can be reached by decreasing the ratio Ej /Ec at 
given fixed length N . However this limit is beyond the 
scope of the present work which is based on the single 
QPS approximation. 

Alternatively, for a given ratio Ej / Ec, we can increase 
the length TV of the system. The finite size effects that 
we have discussed in the previous section for the phase 
difference 6 — it vanish as 1/N. We now show that the 
finite-size corrections due to the asymmetry of the barrier 
vanish more rapidly, namely as 1/N 2 . Consequently, in 
the cross-over range in which N spans from N m i n = 6 to 
N S> 1, we can neglect the difference between v as (for 
S 7r) and v (for 5 = n). It is then justified to use the 
renormalized amplitude v for the symmetric potential, 
Eq. (|4"B")) . when calculating the maximum supercurrent 
for the model Eq. (|T3"|) in the full range of q. 

Let us introduce k = 1 — 5/tt as the natural parameter 



to quantify the asymmetry. For an asymmetric potential, 
we can approximate the hopping energy between the two 
levels as the geometrical average obtained by considering 
the hopping for the left part of the potential and the right 
part of the potential with respect to the maximum ; 36 ! 37 

v^y/viMR- (49) 

The left (right) amplitude Ul { v r) is the tunneling 
amplitude for the symmetric double potential, Eq. (|48|) . 
with the barrier height 2E*j given by the difference be- 
tween the maximum energy and the left (right) mini- 
mum. In a similar way, AO* /2-k is the distance be- 
tween the maximum and the left (right) minimum point. 
Keeping the leading term of the expansion in 1/N, they 
read E}/E.j = 1 - (tt 2 /4A) ± (n 2 /2N)n and A6*/2n = 
1 — 1/A± (tt/N)h where ± is for L and R. Inserting this 
expansion into Eq. (|4T)]) . we finally conclude that finite- 
size corrections associated with the asymmetry parame- 
ter k cancel at order 1/N. When higher-order corrections 
are taken into account by using the full formula Eq. (|49|) 
for the QPS amplitude v = v{8), the differences with the 
amplitude at v = v(tt) are practically unnoticeable. 



V. EFFECTS OF THE CAPACITANCE Co 

We now consider the effect of the ground capacitance 
Co on the QPS amplitude. In view of the discussion of 
the previous section, we will focus on the analysis of the 
case 6 = tt. Then the effective potential appearing in the 
first line of Eq. Q37p is replaced by a renormalized cosine 
potential 
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(50) 



where E* c , E} are defined in Eqs. (|45I46|) and A6* in 
Eq. (|47|) . Within the semiclassical instanton approach^ 
the QPS amplitude reads 



Aexp 



Sif f /H 



(51) 



where S^J is the effective action Eq. (f5T)]) evaluated at 
9 c i(t), the asymptotic path which minimizes the action 
and which connects the two relevant minima in the limit 
ft — > oo, i.e., the instanton solution. The prcfactor A 
is related to the quantum fluctuations around this mini- 
mum pathi2£ 

In contrast with the previous analysis for Co = 0, we 
now have to take into account the effects of the kernel 
G(t) in the action Eq. ((50)) . 



A. Parabolic approximation 

The first step is to find the classical path c i(t). For a 
cosine-potential, this solution is known analytically only 
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when the kernel is zero G(r) = (Co = 0). For the 
general case, to the best of our knowledge, the solution 
is unknown. We use the Villain approximation to solve 
the problem^ We replace the periodic cosine potential 
by a parabolic potential 



dr 
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(52) 



It is worth noting that, for G = (C = 0) and TV > 1, 
the instanton solution for this potential yields the action 



S 



(ci) 



>/h = ir*(Vj/8E c ) 1 ' 2 ,i.e. 
different from the one found for a cosine potential with 
the same amplitude Vj. Thus, in order to recover the 
previous results for Co = 0, it is convenient to set the 
height of the parabolic periodic potential to the value 
Vj = (8/ir 2 ) 2 Ej to take into account the difference be- 
tween the profiles of the two potentials. 

After the calculation (see Appendix iBl for details), the 
action Eq. (f5"2")l with the instanton path reads^S 



the numerical coefficient is 



= 2nVj ( 1 



TV 



dw- 



1 



C* + iC + 



4e 2 Vj/fi 2 
T 



G(u>) 



(53) 

where G(u>) is the continuous limit for the Fourier trans- 
form G(w e ) = G e defined in Eq. (J3HJ). In Fig. M the 

behavior of Sp is shown function of the ring 

(d) 

size TV. For Co = and N S> 1, Sp saturates to 
So = (SEj/Ec) 1 ^ 2 - When the ground capacitance is re- 
stored, we find a logarithmic scaling with TV for TV > A. 
Specifically, we find that v ~ exp[— a log(TV)] = 1/TV Q 
with a = 7r y/Ej/(8E ) (see Appendix 0. 
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Figure 9: For Ej = 8Ec, the classical action Sp C '' Eq. (JS3J) 
scaled with So = (SEj/Ec) 1 ^ 2 as a function of TV for different 
screening lengths A = 3, 10, 50, oo. Notice the cross-over at 
TV > A to the logarithmic scaling ~ log (TV). 



B. The prefactor A 

So far our semiclassical approach was general. We will 
now restrict ourselves to the calculation of the finite-size 
corrections entering the renormalized amplitude v with 
exponential accuracy; this corresponds to the leading de- 
pendence on TV. Correspondingly, we will use an approxi- 
mate expression for the prefactor A. This prefactor is as- 
sociated with the quantum gaussian fluctuations around 
the classical path. Specifically, we neglect the contribu- 
tion of the low energy paths having a mean kinetic energy 
lower than the height of the potential. 

Proceeding as in Sec. IIII D[ to estimate the kinetic en- 
ergy we determine the effective capacitance of the junc- 
tion in the ring but now taking into account the finite 
size effects. For tag — > oo the kernel Eq. (j38]l reduces to 
Gt = AC(huj e ) 2 /(4e 2 ) where 



AC 



C„ 



2 (TV - 1) 



E 

fc=i 



1 + cos 






1— cos 







(54) 



so that the effective capacitance of the junction corre- 
sponds to 



C 



eff 



C* 



\co + AC. 



(55) 




Again, from Eq. ([52"]) , as Vj ~ Ej we see that the thresh- 
old frequency separating the high and low energy regions 
is still uj£ = uj m ax for C > Co- From these observations, 
in order to take into account the contribution of the high 
energy quantum fluctuations, we replace the effective ca- 
pacitance in the prefactor A of Eq. (jig)) (case Co = 0), 
with the effective capacitance Eq. (|55p . The result reads 



(56) 



We now discuss the conditions for this approximation 
to be valid. For uj e < uj max the kernel G in Eq. ([52")) 
couples the dynamics of the winding phase to the one of 
the modes. As we have seen at the end of Sec. IIII Dl at 
TV = oo, for C 3> Co this interaction between the wind- 
ing phase and the modes corresponds to a perturbation 
at low energies for the Gaussian harmonic fluctuations. 
Indeed, at TV = oo we have weak damping for the qual- 
ity factor Q = 2tt(C/Co) 1/2 > 1. Although at finite TV 
the dynamics of the JJ ring does not correspond to a 
real resistance, the ratio C /Co still plays the role of a 
dimensionless coupling between the winding phase and 
the harmonic bath. The effect of this interaction on the 
low energy Gaussian fluctuations can be neglected in the 
prefactor but not in the exponent where, as we have seen, 
they strongly affect the instanton classical path. 
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C. Superconductor-Insulator transition 



By inserting Eq. (|53| and Eq. (f56| into Eq. (foTj) . wc 
obtain the QPS amplitude which we use to calculate the 
maximum supercurrent. Here we discuss some numerical 
results valid in a general range of parameters while in 
Appendix [C] we discuss some analytical results valid for 
very long chains (TV ^> A). 

In Fig. 1101 we plot the factor q = N 2 u{N) j '{2tt 2 ) as a 
function of the ring size TV for the values C = 2Co and 
for different ratios of Ej/Eq- In the inset of Fig. [TU1 wc 
also show the behavior of vjEj. The first striking ob- 
servation is that for large ratios of Ej/Ec = 2(Ej/Eq), 
the parameter q does not increase as TV 2 but, indeed, it 
vanishes upon increasing the length. This is due to the 
logarithmic dependence of the action S e f / that we have 
obtained in Sec. IV Al 

The presence of the logarithmic dependence is related 
to the lowest energy modes (see Appendix ICj) . Indeed, 
below the corresponding threshold uj£ <C w r nin , the wind- 
ing junction feels the discreteness of the spectrum of the 
environment with which it can exchange energy. It is 
equivalent to say that there is no real dissipation at low 
frequency. An instanton solution conserving the initial 
energy still exists with a finite value of the correspond- 

1/2 

ing action. For TV ^> A, we have uj m i n ~ (8EjEq) ' tt/TV. 
Eventually, w m m vanishes as TV — > oo and the action as- 
sociated to the instanton diverges. 

As a consequence, the parameter q behaves as q ~ 
N 2 v ~ TV 2 ~" and it scales cither to zero or to infinity 
for TV = oo. A cross-over is therefore expected at some 
critical ratio Ej /Eq as a = it \J E j / (SEq) (see Appendix 

EJ). 



1e+02 



1e-00 



1e-02 




1e-08 



10000 



Figure 10: The parameter q as a function of TV for C = 2Co 
and different ratios Ej/Ec- Inset: the behavior of the QPS 
amplitude u scaled with Ej. 



for Ej/Eq < 3.24 it vanishes. On the other hand, we 
remark that for finite systems, the current shows a non- 
monotonic behavior: it increases with iV up to a maxi- 
mum after which it decreases. 
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Figure 11: The maximum supercurrent Imax scaled with the 
classical value 7 c i as a function of TV for C = 2Co and different 
ratios Ej/Ec- 

The interplay between finite size effects, which scale as 
v(N) ~ cxp(N c /N), and the low energy modes at finite 
Co, which reduce the QPS amplitude as v ~ l/N a , cause 
that the total QPS amplitude of the ring v r i ng = Nv(N) 
has a weak iV-dependence so that the dependence of the 
maximum supercurrent on TV appears almost flat over a 
large range of N when the critical ratio is approached 
(see Fig. El). 

Due to the non-monotonicity of the current as a func- 
tion of TV, the critical point between the superconducting 
phase and the insulator phase can be better determined 
by plotting the maximal supercurrent as a function of the 
ratio Ej/Ec for different circumferences. The result is 
shown in Fig. [12] for the case C = 2Cq. 
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The behavior of the maximum supercurrent I m ax 
scaled with I c i as a function of TV, Fig. [TT1 is related to 
the behavior of q. In the very long length limit, I m ax/ Ic\ 
increases and saturates to one for Ej/Eq > 3.24 whereas 



Figure 12: Imax scaled with 7 c i as a function of Ej/Ec = 
IE j / Eo for rings of different circumferences. The dashed ver- 
tical line corresponds to the critical value Ej/Eq = 32/ ty 2 of 
Ref. QJ in the limit C> C . 
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We find that the critical value saturates to the value 
Ej/Eq ~ 3.24, as long as C > 2C . This result is in 
agreement with the result obtained by Korshunov^ who 
calculated Ej/Eq = 32 /ir 2 in the regime C > C . Re- 
ducing the ratio C/Cq, we find a small increase of the 
critical value Ej/Eq which saturates to ~ 4.16 for van- 
ishing mutual capacitance C = 0. In this limit, we do 
not recover the critical ratio Ej/Eq = 128/-7T 2 obtained 
by Bradley and DoniachJ^ This is expected as we will 
now explain. 

D. Discussion 

Mathematically, the onc-dimensional quantum Hamil- 
tonian of the JJ chain can be mapped onto a two- 
dimensional classical spin HamiltonianJ^ For instance, 
when C = 0, the effective model corresponds to the 
two-dimensional classical (anisotropic) XY Hamiltonian. 
Hence, according to this mapping, the superconductor- 
insulator phase transition occurring in the case N = oo 
corresponds to the order-disorder phase transition of 
an ensemble of ferromagnetically coupled planar classi- 
cal spins. This is the celebrated Bcrezinskii-Kostcrlitz- 
Thouless transition, marked by a disruption of the or- 
dered ferromagnetic phase due to the appearance of vor- 
ticity: topological defects for which the spin orientation 
changes by 2ir when following a closed path around them 
once. When Ej Ec, a quantum phase-slip event cor- 
responds to a "vortex" created in the ordered phase in 
the plane defined by the imaginary time axis and the axis 
corresponding to the positions of the islands.— 

In the XY-model, to which the Lagrangian of Eq. p^|) 
can be mapped for C = 0, the interaction between vor- 
tices can not be ignored as the transition is determined 
by the dissociation of vortex - anti-vortex pairs. Con- 
sequently, our single instanton approximation takes into 
account only isolated free vortices; it applies only quali- 
tatively in this limit. 

On the other hand we observe that, for the general 
case of (C, Co), the two-dimensional classical spin Hamil- 
tonian associated with the Lagrangian of Eq. (fl6|) has 
a four-body anisotropic interaction between the spins 
ip n (T). Then the term C /Cq represents the ratio between 
the two coupling constants, the first associated with the 
interaction along the imaginary time axis fi. 2 C/4e 2 and 
the second to the interaction along the chain fi, 2 Co/4e 2 . 
In the limit Cq/C < 1, we have a strongly anisotropic 
spin model. 

In this case it is more convenient to use the phase dif- 
ferences {0 n } to express the Lagrangian as ($b = 0) 

^ = E - * ~ SI E Wnfln. , 

71 ^ ' 71,771 

(57) 

where V n = f_ n (dk/2ir)[l — exp(ifcn)]/[l — cos(fc)] (Vq = 
0), and to consider {6 n } as classical spins and the Joseph- 
son term as a local magnetic field. 



Hence, the superconductor/insulator transition is 
better analyzed theoretically as a dissipative phase 
transition 2 ^ in which the local phase-junction 9 n weakly 
interacts with an external mean field (the external en- 
vironment) associated with the other phase-junctions. 
Indeed the condition Cq/C < 1 enables us to consider 
the QPSs as rare events occurring in different junctions. 
Hence we can use a single instanton approximation (in- 
troduced by Korshunov^i) in which we neglect the in- 
teractions between QPSs occurring at different junctions 
and we consider the dynamics of the surrounding phases 
in the harmonic approximation. Finally, the condition 
Ej > Ec enables us to use the semiclassical instanton 
approach to describe the local QPS. 



VI. SUMMARY AND CONCLUSIONS 

In this work we studied the quantum QPS processes in 
ID Josephson junction rings in the strong Josephson cou- 
pling limit Ej 3> Ec , Eq. In contrast with the previous 
work,— we consider J J rings of finite size in a wide range 
of lengths (N > 5) and with ground capacitance Cq. Wc 
calculated the renormalized QPS amplitude v(N) and we 
discussed its consequence for the maximum supercurrent 
Imax flowing through a JJ ring threaded by a magnetic 
flux. 

For the case Cq = 0, wc found an interplay be- 
tween different finite-size effects which gives rise to a 
non-monotonic behavior of the maximum supercurrent 
/I c \ as a function of the ring size. The critical 
length above which finite size corrections are negligible 
isN c ~(Tr 2 +4)y/8Ej/E c . 

When the ground capacitance is restored, Cq > 0, dis- 
persive modes are possible on the ring which are directly 
coupled to the local winding phase-difference. When 
N —> oo we found that the system converges either to 
a superconducting state with I max /I c \ = l or to an insu- 
lating phase I m ax/Ic\ = 0, depending on the ratio Ej/Eq. 
For C > Cq we found as critical ratio (Ej/Eq) = 3.24, 
in agreement with the previous work of Ref. [l5|. 

Although our analysis was mainly developed for a ring 
of identical junctions, it is also relevant for other sys- 
tems, for instance for a weak winding junction coupled 
to a chain of Josephson junctions in which the QPS is 
prevented, i.e. the so-called fluxoniumiii 

We have discussed the validity of the single QPS ap- 
proximation. The regime Ec > Ej associated with the 
regime of strong interaction between the QPSs is beyond 
the scope of this article. A global phase diagram was 
reported some time ago in Refs. [l4lfl5l in the thermody- 
namical limit but it was based on a perturbative renor- 
malization group analysis. A phase-diagram obtained by 
non-perturbative approaches as well as the behavior of 
finite-size systems for arbitrary ranges of the parameters 
Ej,Eq,Ec and C/Cq constitutes a still open theoretical 
issue. 
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Appendix A: Path integral 

In this Appendix we summarize the main steps for the 
calculation of the path integral Eq. 

First we recall that any periodic function ip n on the 
lattice n = 0, . . . , N — 1 can be decomposed as Eq. (fT8|) 
where the set {<f k } are complex numbers, i.e. <p k = <p k + 
itp k . They are related by the condition tfN-k — 'fit as *Pn 
is real. Using this property, we can write 



ifk=o + (~l)'V fe=f +]T (ip k e l ^ n + cc.) 
' fc=i 



N even 

(Al) 

For N even we have k max = (N/2) — 1 whereas for N odd 
we have k max = (N — l)/2. We can express the Euclidean 
Lagrangians only in terms of the independent variables. 
Setting the action 

Lo Ofc) = + \^i k bj\xi , (A2) 

one can demonstrate that the harmonic Lagrangian 
Eq. Q17p (omitting the constant proportional to S m ) in 
terms of the modes k reads 



T (N) ^0.2 . T , s , o 
= Y^=o+ L o(%=f )+2 

N even 



E 

fc=l 



(A3) 

We now consider the actions Eqs. ([Ml), ©, and (J3TJ) 
where we have N — 1 harmonic variables tp n . 

Inserting Eq. (|Al[) with N replaced by N — 1 in 
Eq. pip , one can obtain that only the imaginary parts 
of the modes are linearly coupled to external forces 
(6(t),8(t)) through the position and through the veloc- 
ity, i.e. Eq. (|33|) . For the notation we set 



1k{r) =¥>a(t), 



and we can write 



£ 3 = £a 

A,- 



h 2 c 



-q k t 



Ej 



-q k t 



(A4) 



(A5) 



From Eq. (|A1|) we observe that the phase <p n =o at n = 
depends only on the real part of the mode {ip k }. There- 
fore the first term of Eq. (|30|). (p„ = o, couples only the 
real parts of the different modes k. As real and imagi- 
nary part are decoupled, the relevant term in Eq. (|32l) in 



which we are interested for the calculation of the effective 
action reduces to 



£ 2 ~ £ +2 £ (\iikql + \ 



;VkU 2 k qt 



(A6) 



k=l 

with £ Q = Ej8f n (N - 1)/(2N 2 ). We express the generic 
periodic path of the partition function as 



qh(r)=q hfi + Y J {lh,l^ lT +c. C : 



(A7) 



in which the Fourier (Matsubara) component at u>i = 

27r£//3 is given by q k . t = (1/13) f drexp(- iuer)qk(r). 
Inserting the expression Eq. (|A7|) into the Lagrangian 
C-2 + £3 Eqs. (|A5|) . (|A6[) and integrating over the time, 
we obtain 



cIt (£ 2 + £ 3 )=£ + 2Y, S k ,o + E 5 M 



(A8) 



k=l 



The first term S k ^o contains only the component at zero 
frequency of 9(t) and q k (r), 



c ^ 2 2, P E J n / Am 

S k ,o = ~H k u k q kfl + —^-qkfiOo . (A9) 

The second terms S ki e contain all the nonzero frequency 
components (£ > 0) 

S k ,e=PHk (wf + cj^) \q k j\ 2 +^^-(uj 2 e +ujl) (ql^Oe + cc.) 

(A10) 

Finally, the path integral is evaluated by integration over 
the Fourier components as 



AZ = f[ <fv[q k (r)} e 
fc=i ^ 



i r4 



/ ( fdr(£ 2 +£ 3 ) 



_££o 

e 



n" ax f dq k ,o e ft -pr 
/ I i-r-RK i- 1 



fc=l 



4^/3 ft 



(=1 



k,e""ik,e _ -°fc,< 
— - e » 



(All) 



where the pre-exponential factors are the Jacobians asso- 
ciated with the transformation of the path integral from 
the time-space to the frequency-spaced (iknlkt) 
note respectively the real and imaginary part of q k .i . The 
integral Eq. (jAllj) has the general Gaussian form and its 
evaluation is straightforward. After some algebra, the 
relevant expontial term is 



AZ 
with 



111 



exp \ -—— V Y(u e ) )-A—fl 



m 2 , (A12) 



Y(ut) = 



N 



sin 2 


' 27Tfe ^ 


) 


l-cos^) 


+ 2^1 









(A13) 
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By adding the exponential term of Eqs. (|A12[) . (|A13|) to 
the Lagrangian C\ Eq. (j2T)]) and going back to the time 
representation, we obtain the result shown in Eqs. (|36|) . 
([57]). and 



Appendix B: Instanton solution in periodic 
parabolic potential 

Referring to the action Eq. (f52"j). we set huij = 
4Vje 2 /(C* + Co/2) as a short-hand notation. The most 
general path can be always expressed as 



'(r)=0(- J 8/2) + dr'9(r') 



(Bl) 



For the asymptotic instanton-like solution, we require 
the following boundary conditions at the end-points 9i = 
6{-P/2) and 9 f = 6(fi/2) 



lim Qi 



0, lim Of 

[3->+oo 



A9* 



(B2) 



It is useful to use the Fourier components of the velocity 
as free variables 



'(r) 



+ 00 

E < 

*=-oo 



ijjre^e-^. (B3) 



By assuming that the total energy is conserved, we im- 
pose that the initial velocity is equal to the final one: 
9(—B/2) = 9(f3/2). More specifically, due to the symme- 
try of the potential, we can assume that the velocity is an 
even function of the time 6(t) = 0a+ J2t=i 2 cos(ui£t)9i, 
with 0£ are real numbers. By definition, the average ve- 
locity is 9q = (Of — Si) I '(3. In terms of the variables {9e}, 
the path Eq. (|B1[) reads 



(r) 



% (r + fi/2) 




lOJi 



(B4) 

We recall the definition of the kernel G(r) in the action 
Eq. ([52]) in the Fourier space 



with G-f = G* p 



(B5) 



Using the expressions Eqs. (|B3I) . (IB4I) .([B5 [) in the action 
Eq. (|52[) and taking into account the symmetry respect 
to the time, the time integration is straightforward. The 
general action is thus expressed in term of the vari- 
ables {9t} 



S P ({0i}) 



Vj 
2w» 



Vj , Vj , Gi 



1=1 ' 1 



+ 



Vj 
2 



a + ^ 

1 2 



+ 



pe 
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(B6) 



Then, by the condition of minimization dS p /d8( = 0, 



we find the classical solution 



S(cO 



2^/(3 



[i - (-i)l e i + ^ 



(B7) 



We insert the solution Eq. (|B7[) into the action Eq. (|B6[) 
to obtain the value of the action at the minimum 



Op 



PV.,91 , AVj 



+ 



[['-<- 



1) 



2 



h 2hui HB 

J ' i=i w t t- i +u 2 Gi/(uiVj) 

(B8) 

Finally, we take the limit B oo. We use the bound- 
ary conditions Eq. (|B2[) and the fact that the average 
velocity 9 scales as (9 f - 9i)/B ~ A9*//3. After that, 
the resulting series with respect to Matsubara frequen- 
cies £ = 1 , . . . , oo converges to an integral uii = us and 
the resulting action coincides with Eq. (|53[) of the main 
text by recalling Huij = 4Vje 2 /(C* + C /2) and Eq. {47]) 
for AO*. 



Appendix C: Long J J ring with ground capacitance 

In this Appendix we discuss some general features 
and some analytic limits of the effective action obtained 
within the parabolic approximation Eq. (|53j) . To simplify 
the notation we set Vj = Ej as the results are qualita- 
tively the same for two different coefficients. We focus 
on the long circumferences limit of the rings defined by 
the condition that we can neglect the corrections of or- 
der (I/AT) in E* c , E} and A9* as well as N > max(A, 1). 
In this regime, we replace the sum with respect to the 
modes k in G(uig) = G(u>) Eq. (|38|) with an integral and 
we obtain 



G(w) 



(hwf 

16e 2 



Co F (ui) J1 + 



AC 



- 1 



(Cl) 



where we set 



F(u) = dl + 



1 arctan ( mm \ 1 1 



(C2) 

where uj m ax, ui m in are defined in Eqs. (|22 p .(|25j ) . Note that 
we can not neglect the N— dependence in the function 
F(u>), Eq. (|C2[) . as it is strongly dependent on the mini- 
mum cut-off frequency uj m i„. Using the Eqs. (|Clj) . (|C2[) . 
we can express the action on the classical instanton path 
Eq. (USD as 



2tt—^- I du) 



(C3) 



W 1+2F (^/(l+lc/Ca) 1 '' 2 



h h J 

The two frequencies Lo max , w„ ljn define three ranges for 
the integral on the frequency uj. They are: i) u? max <C 
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lo, ii) u min < oj < o; TOax , and iii) w < w TOin . In 
these range the function F{uj) can be approximated as: 
i) F(w) ~ 1, ii) ~ u max /u), and iii) ~ 

(N/ir)/(l + 4C/C ) 1/2 , to leading order in 1/JV. Cutting 
the integral Eq. (|C3j) in three parts, we use the three ap- 
proximated expression for the function F(u>) to evaluate 
the integration 



L 



+00 
duj - ■ 



f+00 






/ dui ■ ■ ■ 


+ duj ■ ■ ■ 




' Umax 




Jo 



(C4) 



The first integral i) in the high-frequency range as well 
as the third integral iii) in the low-frequency range gives 
a result independent of the ring's size N. On the other 
hand, for the second integral ii) in the intermediate fre- 
quency range we find the important result 



S 



(ct) 



E., 



2tt- 



ttEj 
8Eq 



du 



LJ 2 +LOLO max (l+4C/C ) 1/2 



log(iV) 



(C5) 



As we have reported numerically, we have a logarith- 
mic dependence on N of the classical action yielding 
to a power-law dependence of the QPS amplitude v. 
This leads to a superconductor/insulator phase transi- 
tion when the limit N = 00 is taken. 



Finally, we observe that the function F(w) Eq. (|C2[) 
saturates to a constant for uj <C kWin This low-frequency 
cut-off is important because it makes that the total in- 
tegral convergent. As we have explained in the text, be- 
low this threshold u> -C i^mim the winding junction feels 
the discreteness of the spectrum of the environment with 
which it can exchange energy. That corresponds to say 
that there is no real dissipation at low frequency and 
an instanton solution conserving the initial energy still 
exists. 
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